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ABSTRACT 

Unified schemes of active galactic nuclei (AGN) require an obscuring dusty torus around the central engine. 
The compact sizes (only a few pc) determined in recent high-resolution observations require that the obscuring 
matter be clumpy and located inside the region where the black-hole gravity dominates over the galactic bulge. 
This location is in line with the scenario depicting the torus as the region of the clumpy wind coming off the 
accretion disk in which the clouds are dusty and optically thick. We study here the outflow scenario within the 
framework of hydromagnetic disk winds, incorporating the cloud properties determined from detailed modeling 
of the IR emission from clumpy tori. We find that torus clouds were likely detected in recent water maser 
observations of NGC 3079. In the wind scenario, the AGN main dynamic channel for release of accreted mass 
seems to be switching at low luminosities from torus outflow to radio jets. The torus disappears when the 
bolometric luminosity decreases below ~ 10 42 ergs -1 because the accretion onto the central black hole can 
no longer sustain the required cloud outflow rate. This disappearance seems to have been observed in both 
LINERs and radio galaxies. With further luminosity decrease, suppression of cloud outflow spreads radially 
inward from the disk's dusty, molecular region into its atomic, ionized zone, resulting in disappearance of the 
broad emission line region at lower luminosities, yet to be determined. 

Subject headings: galaxies: active, galaxies: nuclei, galaxies: Seyfert, infrared: galaxies, quasars: general, 



1. INTRODUCTION 

The great diversity of AGN classes has been explained by a 
single unified scheme (e.g., Antonucci 1993; Urry & Padovani 
1995). The nuclear activity is powered by a supermassive 
(~10 6 -10 10 M Q ) black hole (SBH) and its accretion disk. 
This central engine is surrounded by dusty clouds, which are 
individually optically thick, in a toroidal structure (Rrolik & 
Begelman 1988). The clumpy torus is characterized by two 
fundamental properties: (1) Anisotropic obscuration of the 
central region so that sources viewed face-on are recognized 
as type 1 objects, those observed edge-on are type 2. From 
the statistics of Seyfert 1 and 2 galaxies, the torus height 2 
and radius obey H/R ~ 1 (Schmitt et al. 2001). (2) Dust 
re-emission in the IR of the AGN obscured radiation. This 
emission is the only means to find the torus radial extent since 
obscuration cannot determine separately H and R. 

Two approaches have been taken for the torus dynamic ori- 
gin. A hydrostatic scenario depicts the torus as a doughnut- 
like structure populated by molecular clouds accreted from 
the galaxy (Rrolik & Begelman 1988). However, the origin 
of vertical motions capable of sustaining the clouds in a hy- 
drostatic structure with H ~ R was recognized from the start 
as problematic and has eluded solution thus far (e.g., Davies 
et al 2006). The other scenario, based on the seminal work by 
Blandford & Payne (1982), involves the outflow of clouds em- 
bedded in a hydromagnetic disk wind (Emmering et al. 1992, 
EBS hereafter; Konigl & Kartje 1994; Kartje & Konigl 1996; 
Bottorff et al. 1997, 2000; Kartje at al. 1999; Everett 2005). 
In this approach the torus is merely a region in the wind which 
happens to provide the required toroidal obscuration, i.e., it is 
that region wherein the clouds are dusty and optically thick. 

Recent high-resolution IR observations indicate that the 
torus size might be no more than a few pc (Elitzur 2005, and 
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references therein); in particular, VLTI observations of NGC 
1068 show that the FWHM size of the 12 fim. emission is only 
~ 4 pc (Jaffe et al 2004). The compact sizes place the torus 
inside the region where the SBH gravity dominates over the 
galactic bulge. Consider a bulge that induces a linearly rising 
rotation curve with £1 ~ lkms -1 pc -1 , as is typical of AGN 
host galaxies (Sofue et al. 1999), and an SBH with a mass of 
M, 7 x 10 7 Mq at its center. The SBH will dominate the grav- 
itational motions within a radius 35 pc x (M.r/nj) 1 / 3 , where 
ill = fl/(lkms -1 pc -1 ). Since the torus is well within this 
region, its dynamic origin is determined in all likelihood by 
the central engine and its accretion disk, giving strong impetus 
to the outflow paradigm. The central accretion disk appears 
to be fed by a midplane influx of cold, clumpy material from 
the main body of the galaxy (Shlosman et al 1990, and refer- 
ences therein). Approaching the center, conditions for devel- 
oping hydromagnetically- or radiatively-driven winds above 
this equatorial inflow become more favorable. The disk-wind 
geometry provides a natural channel for angular momentum 
outflow from the disk and is found on many spatial scales, 
from protostars to AGN (Blandford & Payne 1982; EBS; Fer- 
reira 2006). Although theoretical questions involving cloud- 
uplift and wind-driving mechanisms remain (e.g., EBS; Fer- 
reira 2006), observations give ample evidence for winds and 
cloud motions (e.g., Arav et al 1997; Elvis 2004) and support 
for a disk-wind geometry (e.g., Hall et al. 2003) in AGN. 

dumpiness provides a natural explanation for the compact- 
ness of the IR emission from AGN tori, and detailed model- 
ing of this emission sets phenomenological constraints on the 
torus cloud properties (Nenkova et al. 2002, 2006; Elitzur et 
al. 2004; Elitzur 2005). Here we explore whether the torus 
outflow scenario can properly accommodate these constraints, 
and the corollaries for the AGN paradigm. 

2. CLOUD PROPERTIES 

The only property of individual clouds constrained from the 
IR modeling is their optical depth. It should lie in the range 



2 



Elitzur & Shlosman 



~ 20-150 at the V band, i.e., column densities iV H ~ 10 - 
10 23 cm -2 assuming standard dust-to-gas ratio. Clouds up- 
lifted into the wind expand while moving away from the disk. 
A cloud starting with 7V H ~ 10 23 cm -2 ceases to partake in 
obscuration when its column is reduced by ~ 100. A cloud 
starting with a smaller column will rise to a smaller height 
before leaving the obscuring torus. Consider the tidal torque 
of the central SBH at a distance r pc in pc, where the Kep- 

lerian period is tx = 3x10 M. 7 r pc yr. To prevent a 
cloud with density p — m p n from shearing, it must be at 
least partially confined by its own gravity and/or the ambi- 
ent magnetic field B. In the former case, the characteris- 
tic Jeans timescale is tj ~ (Gp)^ 1 ^ 2 = 3xl0 4 n 7 1 ^ 2 yr, 
where n-j — n/(10 7 cm -3 ). Therefore, resistance to tidal 
shearing (ij < £k) requires n-j > M, 7 /r pc , which, for a 
cloud with column Nh = A^H,23xl0 23 cm -2 , implies a size 
a < 10 16 A^H,23^p C / M, 7 cm. The corresponding cloud mass 
is M c] < 7xl0" 3 iVH,23a 2 6 M 0) where a w = a/(10 16 cm). 
The self-gravity of such a low mass cannot hold it together 
against dispersal by any reasonable internal velocities. How- 

1/2 

ever, an external magnetic field B ~ 1.5 05 n 7 ' mG would 
suffice if the internal velocity dispersion is 0-5 x 1 km s _1 . 

Clouds with these very same properties may have been de- 
tected by Kondratko et al. (2005) in recent H2O maser obser- 
vations of NGC 3079. Individual maser features detected in 
this source have sizes < 2 x 10 16 cm and span the radial range 
~ 0.4-1.3 pc. Similar to other objects, most maser features 
in this source reside in an edge-on rotating disk. However, 
for the first time, four features were found significantly out 
of the disk plane yet their line-of-sight velocities reflect the 
velocity of the most proximate side of the disk. Kondratko et 
al. note that this can be explained if, as proposed by Kartje et 
al. (1999), maser clouds rise to high latitudes above the rotat- 
ing structure while still maintaining, to some degree, the rota- 
tional velocity imprinted by the parent disk. Because the de- 
tected maser emission involves cloud-cloud amplification that 
requires precise alignment in both position and velocity along 
the line-of-sight, the discovery of four high-latitude maser 
features implies the existence of many more such clouds par- 
taking in the outflow. We suggest that the high-latitude masers 
are yet another manifestation of the dusty, molecular clouds 
that make up the torus region of the disk-wind. Densities of 
H 2 masers are frequently quoted as ~ 10 8 -10 9 cm~ 3 , but 
these are the optimal densities to produce the highest possible 
brightness from a single cloud. Because of the scaling proper- 
ties of H 2 pumping (Elitzur et al. 1989), maser clouds with 
n ~ 10 7 cm~ 3 and TVh ~ 10 23 cm~ 2 produce near-optimal 
inversion and detected radiation in cloud-cloud amplification. 

To produce a torus with the required thickness, the column 
density of a cloud must remain above ~ 10 21 cm~ 2 during 
the rise to height H <~ r. If the cloud expands at the speed 
of sound c s while rising with velocity v, then its size will in- 
crease by a factor x ~ 1 + (r/a)(c s /v). We expect c s /v ~ 
10 ~ 2 since c s is 1 kms -1 at a temperature of 100 K and the 
velocity scale of the cloud outflow is expected to be compara- 
ble to the Keplerian velocity, 208(M. 7 /r pc ) 1/2 kms" 1 . With 
a/r ~ 10~ 3 , the 3-D expansion of such a cloud will cause its 
size to increase by a factor x ~ 10 and the column density 
(oc M c \/a 2 ) to decrease within the torus region by a factor 
x 2 ~ 100, as required. Furthermore, in a cloud-loaded hy- 
dromagnetic disk wind the ambient magnetic field suppresses 
cloud expansion across the field lines while guiding the cloud 



motion. Magnetic confinement can be efficient in either 3-D, 
e.g., the "melon seeds" model of diamagnetic clouds (Rees 
1987), or 2-D, e.g., the slender tubes model of paramagnetic 
clouds (EBS). In the simplest picture of constant area flux 
tubes and thermal expansion only along the field lines, the 
cloud column density remains unchanged in the longitudi- 
nal direction while decreasing as 1/x in the orthogonal di- 
rection. In reality, the decrease of column density is likely 
to fall somewhere between the 1/x 2 of 3-D thermal expan- 
sion and 1/x of constant field. Whichever the case, magnetic 
effects always supplement the thermal confinement, enabling 
smaller and denser torus clouds. 

3. MASS OUTFLOW 

At a bolometric luminosity L = L45 x 10 45 erg s _1 , 

the inner radius of the torus cloud distribution is R[ ~ 

1/2 — ? fi 

0.4pcxL 4 5 T 1500 , where Ti 50 o is the dust sublimation tem- 
perature normalized to 1500 K. This distance was determined 
from the dust temperature on the illuminated face of an opti- 
cally thick cloud of composite dust representing standard mix 
of Galactic grains (Nenkova et al 2006). Some dust compo- 
nents can exist closer than i?;, and attenuation of the central 
UV radiation could also reduce _R; at low latitudes. The torus 
outer radius is YR l7 with Y ~ 5-10. The quantity constrained 
by IR modeling is the number of clouds per unit length, 
-^ci = "ci^ci, where n c \ is the number of clouds per unit vol- 
ume and A c \ the cloud cross-sectional area. Neither n c \ nor 
A c \ enter independently into the calculation of the torus emis- 
sion. An analytic expression for N c \(r, j3) in terms of radial 
distance r and angle j3 from the equatorial plane that agrees 
reasonably with all observations, and can serve as a "stan- 
dard" cloud distribution, is N c \ oc A/"o r~ 2 cxp(— /3 2 /cr 2 ). 
Here A/"o (~ 5-10) is the total number of clouds encountered, 
on average, along a radial equatorial ray, and a (~ 45°) mea- 
sures the angular width of the cloud distribution. This phe- 
nomenological parametrization of the torus cloud distribution 
was not meant to account for the underlying disk, where the 
clouds originate, and should hold only at inclinations j3 above 
the disk surface. In practice, we neglect the disk thickness and 
utilize this expression at all (3. 

The mass density of clouds is M c \n c \= iti p NhN c i, and the 
cloud outflow rate generating the torus is m p Nn J N c \ v c \ dA. 
The integration is over the disk area where the torus clouds are 
injected out of the plane with velocity v c \. This velocity can be 
estimated from the dispersion velocity of molecular material 
in the disk — an outflow is established when the velocity of 
the ordered motion becomes comparable to that of the local 
random motions. From the maser observations in NGC 3079 
Kondratko et al. (2005) find that the velocity dispersion in a 
small region (< 5xl0 16 cm) of strong emission is only ~ 14 
kms -1 , which can be considered typical of the local random 
motions. Therefore, we parameterize the injection velocity as 
v c i{r) = lOkrns -1 xv & u(r) where u = v(r)/v(Ri) and vq is 
the velocity at the inner radius R{ normalized to 10 kms -1 ; 
both v & and u(r) are expected to be of order unity. The torus 
mass outflow rate in the form of clouds is thus 

M w , cl = 0.02 L% 2 JVU^ v 6 x h M yr" 1 . (1) 

Here -/Vj^ = J\foN H ,23 is the column density through all 
clouds in the equatorial plane in units of 10 23 cm~ 2 and 

Ii = J 1 Y (u/y)dyxY/(Y — 1) is an unknown factor of or- 
der unity. The mass inflow rate across the torus is related to 
the bolometric luminosity via M acc = 0.02 L^/e M Q yr -1 , 
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where e is the accretion efficiency at the torus radii. Note that 
M acc is a function of r, since the mass loss in a hydromag- 
netic wind can be a significant fraction of the disk accretion 
rate; at radii characteristic of the torus, M acc can be substan- 
tially larger, and e substantially smaller, than in the innermost 
accretion disk (EBS). The results for M w , c \ and AT acc yield 

Mw.cl Ar tot r -l/2 T _2.6 j n . 

When L decreases, M WiC i/M acc increases. However, this ra- 
tio cannot exceed unity. Moreover, the disk wind applies a 
torque on the underlying accretion disk because each mag- 
netic field line can be considered rigid from its footpoint ro 
and up to the Alfven radius, ta- For a total outflow rate M w 
of a continuous wind loaded with clouds, angular momen- 
tum conservation gives 7 = M w /M acc ~ (fo/^Vi) 2 ' w i tn 
7 ~ 0.1-0.25 either constant (EBS) or weakly dependent on 
r (Pelletier & Pudritz 1992). Therefore, M WiC] < 7M acc . 
When the luminosity is decreasing, eventually this bound is 
violated by eq. |2] and the system can no longer sustain the 
cloud outflow rate required by the torus wind. That is, the 
torus should disappear in low-luminosity AGN. If we take 
7 ~ 0.2, an accretion efficiency e ~ 0.01 (EBS) and all 
other parameters ~ 1, the torus disappearance should oc- 
cur at bolometric luminosities below ~ 10 42 ergs -1 . The 
timescale for this is short. The overall mass in torus clouds 
is Mtot = WpA^H / A cl dV ~ 10 3 N^ 23 L 45 YM Q , so the 

depletion time scale when M w , c i — 7 Macc is M to t/M w ^ c \ ~ 
5xl0 4 (e/j)YN^% 3 years. Since er/7 ~ 1, the torus will 
disappear within a few Keplerian orbits. 

The kinetic luminosity Ek of the torus cloud outflow is gen- 
erally insignificant. A calculation similar to that for M w ,ci 
yields E k = 7xl0 35 L^ 2 N$* 23 v%x I 3 ergs -1 , where I 3 = 

f^(u 3 /y)dyxY/(Y — 1) is of order unity. This is negligible 
in the overall energy budget of the AGN. Finally, the volume 
averaged density of the clouds is <pn, where <fi — aN c \ is their 
volume filling factor. The initial outflow velocity of clouds is 
roughly 10 times larger than their internal velocity dispersion. 
Since we have </> ~ 10~ 2 , the average kinetic energy is similar 
for the cloud outflow and the internal motions. Therefore the 
assumption that the field limiting the clouds expansion is also 
guiding their motion is self-consistent. 

4. SUMMARY AND DISCUSSION 

The analysis presented here shows that the disk wind sce- 
nario is compatible with cloud properties inferred from mod- 
els of the torus IR emission. Furthermore, the very same 
dusty, molecular clouds uplifted from the disk appear to have 
been detected in water maser observations of NGC 3079 
(Kondratko et al. 2005). Proper-motion measurements, if pos- 
sible, of these high-latitude masers would be extremely valu- 
able. Detection of maser polarization seems less likely be- 
cause it would require proper alignment of the magnetic fields 
among maser clouds that amplify each other (Elitzur 1996). 
Co-location of water masers and the torus seems to be ob- 
served also in NGC 1068, where, as noted by Poncelet et al. 
(2006), the two are at similar radial distances. The masers 
there trace a thin, edge-on rotating disk between 0.65 and 1.1 
pc (Greenhill & Gwinn 1977; Gallimore et al. 2001). It is 
intriguing that molecular structures somewhat resembling the 
maser disk, and which might be its extension, are seen also 



on much larger scales in NGC 1068. Galliano et al. (2003) 
find that H2 and CO emissions come from clouds in a disk 
close to edge-on with radius 140 pc and scale height 20 pc, 
for H/ i? ~ 0. 15. From CO mapping, Schinnerer et al. (2000) 
conclude that H ~ 9-10 pc at radial distance of 70 pc, for a 
similar H/R. Detailed, high-resolution mapping of the entire 
region between ~ 1 and 50 pc around the nucleus can shed 
light on the variation in the disk structure during the transi- 
tion from gravitational domination by the galactic bulge to 
the black hole, and on the origin of the disk wind. 

Detectable only in edge-on AGN, the water masers reside in 
the inner molecular regions of the accretion disk. At smaller 
radii, the disk composition switches from dusty and molecu- 
lar to atomic and ionized (e.g., EBS; Kartje et al 1999). From 
time variations of X-ray obscuration that indicate cloud mo- 
tions across the line-of-sight, Risaliti et al. (2002, 2005) find 
that some of the obscuring clouds are at radii considered typ- 
ical of the AGN broad-line region (BLR). Clouds at the very 
inner edge of the molecular disk are protected there from the 
central UV. Upon injection into the wind and exposure to the 
AGN radiation, these clouds can be responsible for the broad 
emission lines at low altitudes above the disk, and for the 
warm absorption at higher altitudes. Indeed, the kinematics 
of such clouds has been shown to fit nicely the BLR line pro- 
files (EBS). Clouds uplifted further out will have their entire 
trajectories at r > i?i and form the obscuring torus. Since 
cloud properties must span some range, the torus structure is 
likely stratified. Clouds that reach the maximal height H ~ r 
must start with a size of ~ 10 16 cm and column density of ~ 
10 23 cm -2 . Clouds with smaller initial size or column will 
lose the required obscuration at lower heights. 

The fraction of type 2 objects appears to decrease with L, in 
accordance with the "receding torus" model (Simpson 2005). 
A decrease of H/R with L arises naturally in the torus out- 
flow scenario. The torus radial dimensions are set from its 
inner radius i?j and increase with L due to dust sublimation. 
On the other hand, since outflow velocities are expected to 
decrease with r, thermal expansion will reduce the cloud col- 
umn density below the range of torus clouds at H < r if the 
relative size a/r is roughly constant. An additional possible 
factor limiting the height at larger L is bending of the stream- 
lines by radiation pressure (e.g., de Kool & Begelman 1995). 

A key prediction of the wind scenario is that the torus disap- 
pears at low bolometric luminosities (< 10 42 erg s" 1 ) because 
mass accretion can no longer sustain the required cloud out- 
flow rate, i.e., the large column densities. Observations seem 
to corroborate this prediction. In an HST study of a complete 
sample of low-luminosity (< 10 ergs -1 ) FR I radio galax- 
ies, Chiaberge et al. (1999) detected the compact core in 85% 
of sources and argued that this high detection rate implies the 
absence of an obscuring torus. This finding was enhanced 
when Whysong & Antonucci (2004 and references therein) 
demonstrated that M87, one of the sample sources, indeed 
does not contain an obscuring torus by placing stringent limits 
on its thermal IR emission; a similar conclusion was reached 
by Perlman et al. (2001). Additional evidence comes from 
UV monitoring of LINERs with L < 10 42 erg s^ 1 by Maoz et 
al. (2005). By detecting variability in most of the monitored 
objects they demonstrated that the AGN makes a significant 
contribution to the UV radiation in each source and that it 
is relatively unobscured in both type 1 and type 2 LINERs. 
The histogram of UV colors shows an overlap between the 
two populations, with the difference between the peaks cor- 
responding to dust obscuration in the type 2 LINERs of only 
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~ 1 magnitude in the R band. Such toroidal obscuration is 
minute in comparison with higher luminosity AGN. 

The predicted torus disappearance at low L does not im- 
ply that the cloud component of the disk wind is abruptly 
extinguished, only that its outflow rate is less than required 
by the "standard" IR emission observed in quasars and high- 
luminosity Seyferts. When M w c i drops below "standard" 
torus values, which at each AGN might occur at somewhat 
different L because of the intrinsic spread of parameters in eq. 
|2| (e, v etc), the outflow still provides toroidal obscuration as 
long as its column exceeds ~ 10 21 cm -2 . Indeed, Maoz et al 
find that some Liners do have obscuration, but much smaller 
than "standard". Line transmission through a low-obscuration 
torus might also explain the low polarizations of broad Ha 
lines observed by Barth et al (1999) in some low luminosity 
systems. The properties of thermal IR emission from such 
low luminosity AGN should differ from the "standard", a pre- 
diction that should be tested in observations. 

If the toroidal obscuration were the only component re- 
moved from the system, all low luminosity AGN would be- 
come type 1 sources. In fact, among the LINERs monitored 
and found to be variable by Maoz et al (2005) there were 
both sources with broad Ha wings (type 1) and those with- 
out (type 2). Since all objects are relatively unobscured, the 
broad line component is truly missing in the type 2 sources in 
this sample. Similarly, Laor (2003) presents arguments that 
some "true" type 2 sources, i.e., having no obscured BLR, do 
exist among AGNs with L < 10 42 ergs -1 (see also Ho et al 
1997). Both findings have a simple explanation if when L 
decreases further, the suppression of mass outflow spreads ra- 



dially inward from the disk's dusty, molecular region into its 
atomic, ionized zone. Then the torus disappearance, i.e., re- 
moval of the toroidal obscuration by the dusty wind, would be 
followed by a diminished outflow from the inner ionized zone 
and disappearance of the BLR at lower, still to be determined, 
luminosities. 

Within the framework of a hydromagnetic wind, M w oc 
M acc , and expressing the luminosity in terms of M acc in eq. 
[l] M w ,ci oc (eMacc) 1 / 2 . The torus cloud and total outflow 
rates both vary in the same direction as the mass accretion 
rate. In contrast, Ho (2002) finds that the radio loudness of 
AGN is inversely correlated with the mass accretion rate (see 
also Greene et al. 2006). That is, when Af acc is decreasing, the 
torus cloud and total outflow rates are decreasing too while the 
radio loudness is increasing. It seems that the AGN switches 
its main dynamic channel for release of accreted mass from 
torus outflow at higher luminosities to radio jets at lower ones, 
with a certain degree of overlap. As noted also by Greene 
et al, X-ray binaries display a similar behavior, switching be- 
tween radio quiet states of high X-ray emission and radio loud 
states with low X-ray emission (see Fender et al. 2004). 
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